Directed flow at midrapidity in ^/snn = 2.76 TeV Pb+Pb coflisions 
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We analyze published data from the ALICE Collaboration in order to obtain the first extraction of 
the recently-proposed rapidity-even directed flow observable v\. An accounting of the correlation due 
to the conservation of transverse momentum restores the factorization seen by ALICE in all other 
Fourier harmonics and thus indicates that the remaining correlation gives a reliable measurement 
of directed flow. We then carry out the flrst viscous hydrodynamic calculation of directed flow, and 
show that it is less sensitive to viscosity than higher harmonics. This allows for a direct extraction 
of the dipole asymmetry of the initial state, providing a strict constraint on the non-equilibrium 
dynamics of the early-time system. A prediction is then made for v\ in Au-Au collisions at RHIC. 

PACS numbers: 25.75.Ld, 24.10.Nz 



I. INTRODUCTION 



Azimuthal correlations between particles emitted in 
heavy-ion collisions are a useful observable to probe the 
behavior of these systems [1] . Specifically, one measures 
the Fourier coefficient [2] 



Ka = (cosnA(^), 



(1) 



where Aiy9 is the relative azimuthal angle between a pair 
of particles, and (• • • ) denotes an average over pairs 
and collisions. The long-range part of this correlation 
(defined by a rapidity gap between the pair) is mostly 
generated by collective, anisotropic flow of the strongly- 
coupled matter created in the collision ^3,. 

The most studied Fourier component is V2A [IHS] 5 cor- 
responding to elliptic flow [7]. Recently it was realized 
that event-by-event fluctuations [2] , generate a whole se- 
ries of harmonics. This has triggered detailed analyses of 
14a forn= 3-6 [BHI^. 

Neglected in these analyses is the first Fourier har- 
monic Via. The observed Via is smaller than V2A and 
V3A [SI in] J and receives a sizable contribution from 
global momentum conservation [TH [TS] which makes its 
interpretation less straightforward. Fluctuations are ex- 
pected to create a dipole asymmetry in the system [TS], 
resulting in a specific directed flow pattern, with high 
transverse momentum particles flowing in the direction 
of the steepest gradient and low pT particles flowing in 
the opposite direction. Hints of this directed flow have 
been extracted from published Via data at the Rela- 
tivistic Heavy-Ion Collider (RHIC) by two of the au- 
thors |17| . and its magnitude and px-dependence were 
shown to be in agreement with ideal hydrodynamic cal- 
culations [TS|. Note that this quantity is distinct from 
the directed flow observable that has been obtained in 
the past from measurements employing a rapidity-odd 
projection |TS]. That rapidity-odd vi gives a negligible 
contribution to Via near mid-rapidity and represents dif- 
ferent physics 20j. 

In this Letter, we show that data on Vla obtained 
by ALICE jTDj can be explained by the superposition of 



two effects: global momentum conservation and directed 
flow. This allows for the first reliable measurement of 
directed flow at mid rapidity. We then carry out the 
first viscous hydrodynamic calculation of directed flow, 
and show that these data can be used to constrain the 
initial dipole asymmetry of the system for each centrality, 
putting strong constraints on models of initial conditions. 



II. DIRECTED FLOW FROM DIHADRON 
CORRELATIONS 

The standard picture of heavy-ion collisions is that an 
approximately thcrmalized fluid is created, which eventu- 
ally breaks up into particles. Particles are emitted inde- 
pendently in each event, with an azimuthal distribution 
that fluctuates from event to event. This yields a two- 
particle correlation which factorizes into the product of 
two single-particle distributions [21^: 



V„a(Pt;-Pt) = Vn[PT)Vn{PT)^ 



(2) 



where the superscripts t and a refer to trigger and asso- 
ciated particles that can be taken from different bins in 
transverse momentum, and w„(pt) is the anisotropic flow 
coefficient. Note that it is possible for the event-averaged 
correlation to not factorize even if independent emission 
holds in each event [21] , and it is also possible for intrinsic 
("non-flow") pair correlations to factorize ,22 . However, 
flow is currently the only known mechanism that pro- 
duces a factorized correlation in the range of transverse 
momentum studied here (the bulk of particles). This 
factorization has been tested in Pb-Pb collisions at the 
Large Hadron CoUider (LHC) [101 [H]: this is done by 
fitting the left-hand side of Eq. (|2|, which is a iV x A^ 
symmetric matrix for N bins in p-p, with the right-hand 
side of Eq. ([2]), using the N values of w„(pt) as fit param- 
eters. The ^LICE collaboration has shown that, while 
the data do factorize for n > 1, this factorization breaks 
down for n — 1 |10j . This is not surprising since there 
is expected to be an additional long-range correlation in- 
duced by momentum conservation that only affects the 
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TABLE I. From left to right: x Psr degree of freedom of 
the fit to the ALICE Via [13 (restricted to pr < 4 GeV/c) 
using Eq. H, and using Eq. ([3|; value of k from the fit; 
estimated value of k from momentum conservation in units of 
10-^(GeV/c)-^ 



first harmonic [T3]. The constraint that all transverse 
momenta add up to yields a back-to-back correlation 
between pairs, which increases linearly with the trans- 
verse momenta of both particles. This correlation adds 
to the correlation from flow; 



Via{PtiPt) = MPt)MPt) - kpTPT- 



(3) 



(Note that the nonflow correlation also factorizes in this 
particular case |52], but the sum does not.) Table ll] com- 
pares the quahty of the fit to Fia using Eq. ([2| or ^. 
Adding one single fit parameter k tremendously increases 
the quality of the fit for all centrality windows. We have 
checked that the values of the fit parameters depend lit- 
tle on the pt window. However, the quality of the fit de- 
creases as higher px particles are included, as observed 
for other harmonics [TD]. Nevertheless, we include the 
entire range of values as a systematic uncertainty in Ta- 
ble|l] varying the lower pT cutoff between 0.25-0.75 GeV, 
and the upper cutoff between 2.5-15.0 GeV. Similarly, we 
use this procedure to estimate a systematic uncertainty 
in vi (see Fig. [2] below) . 

Next, we check whether the value of k from the fit 
is compatible with the value expected from momentum 
conservation. Assuming for simplicity that momentum 
conservation is the only source of correlation, one ob- 
tains [14 k — (X]Pt)~^j where the sum runs over all par- 
ticles emitted in one event, and angular brackets denote 
an average over events in the centrality class. Since ex- 
periments measure only charged particles in a restricted 
phase-space window, only a rough estimate of this quan- 
tity can be made, by extrapolating from existing data. 
We have used the preliminary identified particle px spec- 
tra from ALICE at midrapidity |23j, and extrapolated 
them outside the pT acceptance of the detector using 
Levy fits [24j. In order to extrapolate to all rapidities, we 
have assumed for simplicity that px spectra are indepen- 
dent of rapidity, and we have used the total charged mul- 
tiplicity estimated by the ALICE collaboration ^25j . Neu- 
tral particles were taken into account assuming isospin 
symmetry, and the contribution of particles heavier than 
nucleons was neglected. 

The resulting estimate is shown in the last column of 
Table lU The fit result in general has the correct size and 
increases with % centrality, as expected. The centrality 



dependence is steeper than expected from our rough esti- 
mate, however — the fit value is larger than the estimated 
value for the most peripheral bin, while it is significantly 
smaller for central collisions. We cannot explain this, but 
overall the agreement is reasonable, and a discrepancy in 
k of this size has a very small effect on the extracted 
directed flow; the extracted directed flow curves with k 
fixed to the estimated values were also included in the 
systematic error band in vi, but only have a small effect 
on the two most central bins. 

It has been suggested that the correlation from mo- 
mentum conservation could be larger than our estimate 
because of approximate conservation of transverse mo- 
mentum within smaller subsystems of the entire colli- 
sion system — specifically rapidity slices of roughly unit 
extent [HI [27]. However, we see no evidence here for 
such an enhancement. 

Thus, by taking into account the only obvious non-flow 
correlation, the factorization seen in higher harmonics 
is restored, and we can take the resulting vi{pt) as a 
reliable measurement of directed flow vi (presented in 
Fig. |2] below). 



III. RESULTS OF HYDRODYNAMIC 
CALCULATIONS 

Relativistic viscous hydrodynamics has been shown to 
successfully reproduce u„ for n — 2,3,4 [28.. Here, we 
present the first viscous hydrodynamic calculation for di- 
rected fiow, vi. In hydrodynamics, vi and the corre- 
sponding event-plane angle ^i are defined by vie^^^ = 
(e"^), where angular brackets denote an average over the 
momentum distribution at freeze-out j29| . A collision of 
identical nuclei at mid-rapidity has ip ^ Lp + ir symmetry 
except for fluctuations, hence vi at midrapidity is solely 
due to event-by-event fluctuations in the initial state. 

In event-by-event ideal hydrodynamic calculations, vi 
was found 18 to be approximately proportional to the 
dipole asymmetry of the system ei defined as |16j 



£1 






(4) 



where {•••} denotes an average value over the initial 
energy density after recentering the coordinate system 
{{re^'t'} = 0). 

Here, in order to make a systematic study, we use a 
smooth, symmetric density profile which we deform to 
introduce a dipole asymmetry of the desired size and ori- 
entation. Specifically, our calculation is a 2-1-1 dimen- 
sional viscous hydrodynamic calculation which uses as 
initial condition the transverse energy density (e(r, <j))) 
profile from an optical Glauber model [3D], which is de- 
formed in a way analogous to the previous study of v^, 
and higher harmonics in Ref. |31) : 



j(r, 4>) ^ € frVl + <5cos(0-$i), 0J 



(5) 
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FIG. 1. (Color online) Directed flow wi, scaled by the initial 
dipole asymmetry ei, in a central Pb-Pb collision at 2.76 TeV, 
for different values of the shear viscosity to entropy ratio r\l s. 



where (5 is a small parameter. Both v\ and £\ are pro- 
portional to b for (5 <^ 1. For noncentral collisions, v\ 
depends mildly on the orientation of the dipole asymme- 
try $1 with respect to the impact parameter. Our results 
are averaged over $1. 

Fig. IT] presents the ratio wi/ei as a function of the 
transverse momentum pr for central collisions. Unlike 
higher-order harmonics, which are usually positive for all 
Pt, ^'1 changes sign. The reason is that the net transverse 
momentum of the system is zero by construction, which 
implies (ptVi(pt)) = 0: low-pj^ particles tend to flow in 
the direction opposite to high-py particles. 

The harmonics w„ tend to probe smaller length scales 
with increasing n, and as a result are expected to have an 
increasing sensitivity to viscosity. Our results show that, 
indeed, v\ is less sensitive to viscosity than vi I 32j and 
higher harmonics _28j J51_ . This insensitivity to viscosity 
combined with the approximate proportionality v\ ex e\ 
provides a unique opportunity to place a direct constraint 
on the dipole asymmetry of the early-time system. 

In a realistic Pb-Pb collision, e\ varies from event to 
event. The contribution of directed flow to V\c^ scales 
like e\. Therefore the experimentally measured v\ scales 
like the root-mean-square (rms) value of £\ in the cen- 
trality bin. As we shall see below, there is a wide range 
of predictions for this quantity. With these new data we 
can now quickly discern which are compatible with ex- 
periment by identifying an allowed range of values for the 
rms dipole asymmetry. 

Fig.[2]displays v\ versus pr extracted from ALICE cor- 
relation data using Eq. ([3]). The magnitude and -pT de- 
pendence of v\ are similar at LHC and at RHIC [T7] , and 
the mild centrality dependence, reminiscent of ?;3 [3 [9], 
is expected since both are generated purely from fluctu- 
ations in the initial state. 

The "Pt dependence of v\ in LHC data bears a striking 
resemblance to that predicted by hydrodynamics. Fig. [T] 
In a given centrality window and for a given value of the 
viscosity, one can tune the value of the dipole asymme- 
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FIG. 2. (Golor online) wi (pr) in Pb-Pb collisions at 2.76 TeV 
extracted from correlation data |10| . in various centrality win- 
dows. The shaded band represents the systematic uncertainty 
from the choice of pt window used for the fit. The curves are 
hydrodynamic calculations where the value of ei has been ad- 
justed so as to match the data from above or below, and which 
were used to obtain the upper and lower bound, respectively, 
in Fig. d 



try £1 in the hydrodynamic calculation so as to obtain 
reasonable agreement with data. If one chooses to match 
data at the lowest pr, calculation overpredicts data at 
high pt- Conversely, if one matches data at high pT^ cal- 
culation underpredicts data at low pq^. The correspond- 
ing values of £1 can be considered upper and lower bounds 
on the actual value. 

The values of r\j s (the ratio of shear viscosity to en- 
tropy density) implied by comparisons of elliptic flow 
data to hydrodynamic calculations all lie in the range 
< ry/s < 0.24 [301 |33l E]. Assuming that 77/5 lies in 
this range, we can extract an allowed range for the dipole 
asymmetry, using the extremal values of 77/s and £1 that 
still give a reasonable fit to data, (the extremal curves 
used are shown in Fig. [2]). 

Fig. [3] displays the allowed values of £\ as a function of 
centrality, together with the rms £1 from various Monte- 
Carlo models of initial conditions. The allowed range 
assuming 0.08 < yy/s < 0.16, representing the most com- 
mon values extracted from vi data, are also shown in 
a darker band, to illustrate the small effect of viscos- 
ity. Both the order of magnitude and the centrality de- 
pendence of £1 from Monte-Carlo models resemble the 
allowed values from LHC data. However, there are sig- 
nificant differences between the models. LHC data al- 
ready exclude DIPSY [35] above 10% centrality, and the 
Phobos Glauber model [3S| as well as a recent improved 
mckt model with KNO fluctuations fSSl over the entire 
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FIG. 3. (Color online) Variation of ei with centrality. The 
shaded bands indicate the allowed regions using ALICE data 
in combination with viscous hydrodynamics, assuming either 
< ri/s < 0.24 (lighter shade) or 0.08 < rj/s < 0.16 (darker 
shade). Also shown are the prediction of several Monte-Carlo 
models of initial conditions: DIPSY [35], MC-KLN (obtained 
from mckt vl.OO) (SHOT], mckt vl.25 [371 gS], and two differ- 
ent implementations of the Glauber model utilizing point-like 
nucleons [3H] or uniform disks [3Zj- The binary collision frac- 
tion for both Glauber models was taken to be a; = 0.18 1401. 



Models such as those presented here do not in general 
predict a significant change in dipole asymmetry with 
collision energy. By extracting a best value of ei from 
these LHC data combined with hydrodynamic calcula- 
tions of lower energy collisions, we can make predictions 
for Au-Au collisions at RHIC assuming little change in 
the average dipole asymmetry in a centrality bin. Since 
the change of ei with collision energy predicted by each 
current Monte Carlo model is much smaller than the 
range spanned by the various models, this prediction is 
more reliable than any obtained by assuming a particu- 
lar model for the initial conditions. These are presented 
in Fig. |4J The value of ei at LHC is obtained by taking 
the best fit to the experimental vi for pT < 1.5 GeV/c, 
which is the range where hydrodynamics agrees best with 
data ^9J. Our calculations use rj/s = 0.16 both at LHC 
and at RHIC, but the extrapolation from LHC to RHIC 
depends very weakly on the assumed value of rj/s }30j. 

These predictions are compatible with the attempted 
extraction of vi [17^ from a much more limited set of cor- 
relation data at 20-60% centrality released by the STAR 
collaboration, but a dedicated analysis by one of the ex- 
perimental collaborations at RHIC will allow for a much 
more precise test, including centrality dependence. 



O.lSr 



0.10 



0.05 



0.00 



00-10% 
10-20% 
20-30% 
30-40% 
40-50% 




0.0 0.5 1.0 1.5 2.0 

Pr[GeV/c] 



2.5 



3.0 



FIG. 4. (Color online) Viscous hydrodynamic prediction 
for vi in Au-Au collisions at 200 GeV in various centrality 
windows. 



centrality range. 

Note, however, that since hydrodynamics is expected 
to be more reliable at low transverse momentum, the 
actual value of £i is most likely to lie very close to our 
upper bound. Thus it is possible that the two models 
with the largest dipole asymmetry may only need a slight 
tuning to achieve a correct value, while the value from 
the others may in fact be too low. We prefer here to be 
conservative in our claimed region of allowed values and 
leave it to future study for more stringent conclusions. 



IV. CONCLUSIONS 

We have shown that the first Fourier component of 
the two-particle azimuthal correlation measured at LHC, 
Via, can be explained by collective flow, much in the 
same way as higher harmonics, after the correlation from 
momentum conservation is accounted for. We have thus 
obtained the first measurement of directed flow, vi, at 
midrapidity at the LHC. This experimental result was 
compared with the first viscous hydrodynamic calcula- 
tion of directed flow, vi was found to have a weaker 
dependence on viscosity than V2 and 1)3, which allows for 
the first time a tight constraint to be placed directly on 
the geometry and fluctuations of the early-time system, 
and which rules out certain current theoretical models. 
The extracted values of the dipole asymmetry of the ini- 
tial conditions then allow for predictions to be made for 
directed flow at midrapidity in lower-energy collisions at 
RHIC, which were presented. 
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